RADIATIVE TRANSFER LIMITS OF TWO-FREQUENCY WIGNER 
DISTRIBUTION FOR RANDOM PARABOLIC WAVES 



ALBERT C. FANNJIANG 

Abstract. The present note establishes the self-averaging, radiative transfer limit for the two- 
frequency Wigner distribution for classical waves in random media. Depending on the ratio of 
the wavelength to the correlation length the limiting equation is either a Boltzmann-like integral 
equation or a Fokker-Planck-like differential equation in the phase space. The limiting equation 
is used to estimate three physical parameters: the spatial spread, the coherence length and the 
coherence bandwidth. In the longitudinal case, the Fokker-Planck-like equation can be solved 
exactly. 

Dans cette note nous etablissons la limite auto-moyennisante dans le regime du transfert ra- 
diatif pour la distribution de Wigner a deux frequences dans le cas classique d'ondes en milieux 
aleatoires. Suivant le rapport de la longueur d'onde a la longueur de correlation I'equation limite est 
soit une equation integrale de type Boltzmann soit une equation different ielle dans I'espace d'etat 
du type Fokker- Planck. L'equation limite est utilisee pour estimer trois parametres physiques: 
I'etalement spatial, la longueur de coherence et la largeur de bande coherente. Dans le cas longitu- 
dinal l'equation de type Fokker-Planck admet une solution exacte. 



1. Introduction 

High-data-rate communication systems at millimeter and optical frequencies, remote sensing and 
detection and the astronomical imaging all require understanding of stochastic pulse propagation. 
As pulses consist of a typically broad frequency band the complete information about transient 
propagation requires a solution for the statistical moments of the wave field at different frequencies 
and locations [5:. 

Let ki,k2 be two (relative) wavenumbers nondimensionlized by the central wavenumber /cq- Let 
the wave fields = 1,2, of kj,j = 1,2, satisfy the paraxial wave equation in the dimensionless 
form |S] 

(1) + + ^^^i' = 0' ^' = 

where 7 is the Fresnel number, V represents the refractive index fluctuation with the correlation 
lengths 6 and e in the longitudinal and transverse direction, respectively, and fi is the magnitude. 
Both 5 and e are small parameters related by the anisotropy parameter a as 6a = e. When 
a <C 1 (resp. a ^ 1), the refractive index fluctuates much faster (resp. slower) in the transverse 
direction(s) than in the longitudinal direction. 

An important regime for classical wave propagation takes place when the correlation length is 
much smaller than the propagation distance but is comparable or much larger than the central 
wavelength which is proportional to the Presnel number. This is the radiative transfer regime for 
monochromatic waves described by the following scaling limit 

(2) 7 = fj. = V6, 9 > 0, such that lim 6 < 00, 
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with suitably chosen /i (see jl], jH], 0, ^Uj and references therein). With two different frequencies, 
the most interesting scaHng Umit requires another simultaneous limit 

(3) lim ki = lim k-2 = k, 



lim ko 



Iim-f~^{k2-ki)=f3>0. 



We shall refer to the conditions ((2) and Q as the two- frequency radiative transfer scaling limit. 
The two-frequency mutual coherence function jH] 

Tuiz, X, y) = E[*i(z, X + ^)^2{z, x - ^)], 

where E stands for the ensemble averaging, plays an important role in analyzing propagation of 
random pulses [Hj . But in the radiative transfer regime the two- frequency mutual coherence function 
is not as convenient as the two-frequency Wigner distribution, introduced in JS^, which is a natural 
extension of the standard Wigner distribution and is self-averaging in the radiative transfer regime. 

1.1. Two- frequency Wigner distribution. The two- frequency Wigner distribution is defined 

as 



(4) 



W{z,^,p) 



1 
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ki 2yk\ \/k2 2\/k2 

where, as well as below, we have adopted the easier notation = (kj)^^^^. The choice of the 
scaling factors in the definition is crucial. 

The Wigner distribution has the following properties: 

7y N 
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and hence contains all the information in the two-point two-frequency function. Furthermore, the 
two-frequency Wigner distribution W satisfies the two-frequency Wigner-Moyal equation exactly 

(7) 



dW 



dz 

where the operator Vz is given as 



+ p-V^W 



1 



VM = 



J 2\/ki 



2. Assumptions on the refractive index fluctuation 



We assume that 14 (x) = V{z,x.) is a centered, z-stationary, x-homogeneous random field admit- 
ting the spectral representation 



(x)^/ 



exp (ip • x)Vz{dp) 



with the z-stationary spectral measure Vz(-) satisfying 

nVz{dp)Vz{dq)] = 5(p + q)$o(p)dpdq. 
The transverse power spectrum density is related to the full power spectrum density $(^,p) as 
^o(p) = / ^{w, p)dw. The power spectral density <I>(k) satisfies $(k) = <I>(-k), Vk = {w, p) G 
because the electric susceptibility field is assumed to be real-valued. Hence ^{w,p) = <!>(— t(;,p) = 
— p) = ^{—w, — p) which is related to the detailed balance of the limiting scattering operators 
described below. 
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More specifically we make the following two assumptions. 

Assumption 1. y(z,x) is a Gaussian process with a spectral density ^'(k),k = {w,p) G M'^"'"^ 
which is uniformly bounded and decays at |k| = oo with sufficiently high power of |k|~^. 

We note that the assumption of Gaussianity is not essential and is made here to simplify the 
presentation. 

Let Tz and be the sigma-algebras generated by {Vs : Vs < 2} and {Vg : Vs > z}, respectively 
and let L^iJ^z) and L'^iJ^^) denote the square-integrable functions measurable w.r.t. to them 
respectively. The maximal correlation coefficient r{t) is given by 

(8) r(t) = sup sup E [hg] . 

E[h]=0,E[h2] = l j.[g]^o,E[92] = i 

Assumption 2. The maximal correlation coefficient r{t) is integrable: r{s)ds < 00. 

3. Main theorems 

Theorem 1. Let 6 > be fixed. Then under the two-frequency radiative transfer scaling the 
weak solutions of the Wigner-Moyal equation ^ converges in law in the space C{[0,oo), L'^{M.'^)) 
to that of the following deterministic equation 

(9) ^W + p-VW = [ K{p, q) [e-/5^'i-/(2'=''') W^(x, p + ^) - W{^, p) 
oz 0'^ J I s/k 



where the kernel K is given by 

(10) K(p,q) = ^$(«-l(p + i^).q,q). 

Here and below L^(M^'^) is the space of complex- valued square integrable functions on the phase 
space M^"^ endowed with the weak topology. 

Remark 1. If we now let a — > 0, then the kernel becomes 



i^(P,q) = 6{{p+0=)-q) J <^>{w,q)dw. 



We refer to this as the transverse case because the transverse correlation length e is much shorter 
than the longitudinal correlation length 5. 

On the other hand, in the longitudinal case a — > cxd the limiting kernel would vanish. In order to 
maintain an interesting limit, we increase ^ by a factor of ^/a. Then the kernel for the longitudinal 
case becomes 

(11) i^(p,q) = cl>(0,q). 

In both the longitudinal and transverse cases the fluctuations in the refractive index are extremely 
anisotropic. 

Theorem 2. Assume lime_»o ^ = 0. Then under the two-frequency radiative transfer scaling (0)-|3) 
the weak solutions of the Wigner-Moyal equation converges in law in the space C([0, cxd), L^(]R^)) 
to that of the following deterministic equation 

(12) + p.VW = -k (^iVp - ^x^ • D • (^iVp - ^x^ W{^, p) 
where the diffusion coefficient D is given by 

(13) D(p) = - / ^^,ci)ci^qdq,. 

a J a 
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Remark 2. In the transverse case a ^ 0, the limiting coefficient is 

(14) D(p) = 7r|p|"M <^{w,p±)dw p±(^p±dp±. 

J p px=0 J 

For the longitudinal case a — > oo, we increase fj, by a factor of ^fa as before such that the limiting 
coefficient is nontrivial 

(15) D =7r/$(0,q)q®qdq. 

When ki = k2 or P = 0, eq. @ and (fT^ reduce to the standard radiative transfer equations 
derived in [SJ. The proof of these results fohows exactly the strategy developed in |^ and 
outlined in [0] , originally developed for the standard one- frequency Wigner distribution (see [2] for 
the same strategy applied to the two-frequency Wigner distribution for a different scaling limit). 

Another notable fact is that eq. © with (|1H) and eq. (|12)) with (|15)) are similar to the govern- 
ing equations for the ensemble-averaged two-frequency Wigner distribution for the z-white-noise 
potential analyzed in [S]. This can be understood by the similar behaviors of the potential more 
rapidly fluctuating in z to the z-white-noise potential and the (less) rapid fluctuation in x gives 
rise to self- averaging which is lacking in the z-white-noise potential. 

4. The longitudinal and transverse case 

To illustrate the utility of these equations, we proceed to discuss the two special cases in three 
dimensions. For simplicity, we will assume the isotropy of the medium in the transverse coordi- 
nates such that ^{w,p) = ^{w, |p|). As a consequence D = DI with a constant scalar D in the 
longitudinal case whereas in the transverse case D(p) = C|p|~^p_L (8)p_L with the constant C given 
by 



C = ^jj '^{w,p±)dw\p±\'^dpj 



Here p_L G M? is an unit vector normal to p G M^. 

First of all, the equation (|12|) by itself gives qualitative information about three important pa- 
rameters of the stochastic channel: the spatial spread 0"=,,, the coherence length ic and the coherence 
bandwidth /3c, through the following scaling argument. One seeks the change of variables 

(16) ^ = ^' P = P^cVk, ^ = 7> P = 4- 

where L is the propagation distance to remove all the physical parameters from H12() and to aim 
for the form 

(17) l-W + p-Vs,W = - 
in the longitudinal case and the form 

(18) l^w + p-Vs.W = - f -zVp + ^x) . . ( -zVp + h\W 




d 



z 




in the transverse case. From the left side of (|12j) it immediately follows the first duality relation 
lc<y* ~ L/k. The balance of terms inside each pair of parentheses leads to the second duality 
relation (dc ~ d-Jt/cy^. Finally the removal of or C determines the spatial spread a^, which has a 
different expression in the longitudinal and transverse case. In the longitudinal case, ~ 
whereas in the transverse case a^, ~ 

Cl/3^4/3^-l/6_ 
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In the longitudinal case, the inverse Fourier transform in p renders eq. (|17() to the form 




which can be solved exactly and whose Green function is [3] 

(1 + i\y-y'\'' . (y-y')-(x-x') 



./|2 



(27r)rf sin"'/^ [^1/2(1 



i^cot{V!s(i+i))|y-/3x/2 



y'-;3xV2 



y'-/3x'/2|^an{VI(l+i)) 



xe2 



cos 



for z = 1. This solution gives asymptotically precise information about the cross-frequency corre- 
lation, important for analyzing the information transfer and time reversal with broadband signals 
in the channel described by the random Schrodinger equation jTj (see also jS], CP). It is unclear 
if the transverse case is exactly solvable or not. 
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